#Exercicios 5
##Francisco M. Ulloa S.

esaligna <- read.csv("esaligna.csv", header=TRUE, as.is=FALSE)

#5.1 Editando alguns parmetros grficos

##Crie um grfico de disperso entre dap e ht com:

plot(esaligna$dap~esaligna$ht)

#5.1.1 Legendas dos eixos com nomes das variveis e suas unidades

str(esaligna)
head(esaligna)

plot(esaligna$dap~esaligna$ht, xlab= "Dimetro  Altura do Peito(DAP) (cm)", ylab="Altura Total(m)")


#5.1.2 Marcaes do eixos (ticks) para dentro da rea do grfico

plot(esaligna$dap~esaligna$ht, xlab= "Dimetro  Altura do Peito(DAP) (cm)", ylab="Altura Total(m)", tcl=0.5)


#5.1.3 Apenas dois eixos (formato L)

plot(esaligna$dap~esaligna$ht, xlab= "Dimetro  Altura do Peito(DAP) (cm)", ylab="Altura Total(m)", bty="l", tcl=0.5)

#5.1.4 Ttulo informativo

plot(esaligna$dap~esaligna$ht, xlab= "Dimetro  Altura do Peito(DAP) (cm)", ylab="Altura Total(m)", bty="l", tcl=0.5, main="Grafico de dispersao do DAP(cm) vs Altura Total (m)")

#5.1.5 Tamanho das fontes maiores que o padro

par(cex=1.5)
plot(esaligna$dap~esaligna$ht, xlab= "Dimetro  Altura do Peito(DAP) (cm)", ylab="Altura Total(m)", bty="l", tcl=0.5, main="Grafico de dispersao do DAP(cm) vs Altura Total (m)", cex=1)


####################################################################################################

#5.2 Dois grficos juntos

#5.2.1 Use as variveis dap e talhao para construir dois grficos, colocando-os lado a lado. 

unique(esaligna$talhao)
par(mfrow=c(1,2),mar=c(4, 4, 1.5, 1.5))

##O primeiro deve ser um grfico de desenho de caixa (boxplot) da varivel dap em funo do fator talho. 
## segundo deve ter apenas a mdia e uma barra de desvio-padro do dap, para cada talho.
#5.2.2 Insira tambm uma letra para dizer qual  o grfico a e qual  o b (tanto faz, quem  um e quem  outro).

boxplot(esaligna$dap~esaligna$talhao, ylab="DAP",xlab="Talhao", col = "antiquewhite3")
mtext("a",side=2)

mean(esaligna$talhao)

media <-aggregate(esaligna$dap,by=list(esaligna$talhao),FUN=mean)
sd <-aggregate(esaligna$dap,by=list(esaligna$talhao),FUN=sd)

media$Group.1<-as.factor(media$Group.1)
#mudou a varivel talhao em va. categrica, para poder fazer o boxplot

media_sd<-data.frame(media,sd[,2])
media_sd
colnames(media_sd)<-c("Talhao","Media","Desvio.Padrao")

plot(media_sd$Media,ylim=c(5,21),xlab="Talho",ylab="Mdia do DAP",xaxt="n", pch=8)
axis(side=1,at=1:6,labels=c("16","17","18","22","23","32"))
arrows(x0=c(1:6), y0=(media_sd$Media-media_sd$Desvio.Padrao), y1 = (media_sd$Media+media_sd$Desvio.Padrao),angle=90,code=3)
mtext("b",side=2)

#################################################################################################


5.3 Adivinhando o cdigo

dados53 <- read.csv("5.3.csv", header=TRUE, as.is=TRUE)
str(dados53)

par(mfrow=c(1,2))

str(dados53)
dados53$y2=as.character(dados53$y2)
str(dados53)

plot(dados53$x1,dados53$y1,xlab="",las=1,ylab="",pch=17,xlim=c(0.5,2.2),ylim=c(0,3),bty="l",tcl=0.3,family="serif")
text(x=2.2,y=3,"a",cex = 1.5)

mtext("Euclidean distances",2, cex=1.3, line=2)
mtext("Log(Patch size)(ha)",1, cex=0.8, line=2)
segments(x0=0.5, y0=2.3, x1 =2.1, y1 =1 ,lwd=4)

boxplot(dados53$x2~dados53$y2,xaxt="n",ylim=c(0,3),xlim=c(0,7), bty="l", tcl=0.5)
axis(side=1,at=1:6,c("Small","Medium \n Edge","Medium \n Interior","Large \n Edge","Large \n Interior","Control"))

par(bty="l", tcl=0.3, family="serif")
text(c(1:2,4),y=3,"*",cex=1)
text(x=3,y=3,"**",cex=1)
text(x=5,y=3,"***",cex=1)

tapply(dados53$x2,dados53$y2,max)
which(dados53$x2==1.7208396)
text(x = 6.5, y = 3, "b", cex = 1.5)
dados53[27,3]=NA
